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We investigate the (partial) chiral restoration at finite temperature (T) under the strong external 
magnetic field B — Boz of the SU(2) light-flavor QCD matter. To this end, we employ the instanton- 
liquid QCD vacuum configuration accompanied with the linear Schwinger method for inducing the 
magnetic field. The Harrington-Shepard caloron solution is used to modify the instanton parameters, 
i.e. the average instanton size (p) and inter-instanton distance (R), as functions of T. In addition, 
we include the meson-loop corrections (MLC) as the large- N c corrections because they are critical 
for reproducing the universal chiral restoration pattern. We present the numerical results for the 
constituent-quark mass as well as chiral condensate which signal the spontaneous breakdown of 
chiral-symmetry (SBxS), as functions of T and Bo- From our results we observe that the strengths 
of those chiral order parameters are enhanced with respect to Bo due to the magnetic catalysis effect. 
We also find that there appears a region where the u and d-quark constituent masses coincide with 
each other at eBo « (7 ~ 9)m^, even in the presence of the explicit isospin breaking (m u ^ vm). 
The critical T for the chiral restoration T c tends to shift to the higher temperature in the presence 
of the Bo for the chiral limit but keeps almost stationary for the physical quark mass case. The 
strength of the isospin breaking between the quark condensates is also explored in detail by defining 
the ratio 1Z = ((iu^u) — (id!* d)) / ({iv? u) + (iSd)), which indicates the competition between the 
explicitly isospin breaking effect and magnetic catalysis effect. We also compute the pion weak- 
decay constant F w and pion mass m n below T c , varying the strength of the magnetic field, showing 
correct partial chiral restoration behaviors. Besides we find that the changes for the F n and due 
to the magnetic field is relatively small, in comparison to those caused by the finite T effect. 
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I. INTRODUCTION 

Studies on the breakdown of symmetries and their restorations have been very useful in the analysis of phenomena 
related to phase transitions. Such studies applied to Quantum Chromodynamics (QCD) are particularly fascinating 
since QCD owns a very complicated phase structure. Among them, the (partial) restoration of chiral symmetry 
at finite temperature (T) and/or quark chemical density (p) has been one of the most interesting and stimulating 
subjects for decades. A great number of theoretical and experimental endeavors have been devoted to this subject. 
In particular, the high-T and low-/x region, which resembles the early universe, shows a very noble feature of the 
restoration pattern. Namely there is a crossover phase transition for the nonzero current-quark mass (to/) and that of 
the second-order for to/ = 0. These distinctive patterns of the phase transition are consistent with the universal-class 
argument of the three-dimensional Ising model [1], and turns out to be highly nontrivial in QCD [2, 3]. Moreover, the 
critical endpoint (CEP) [4, 5], tri-critical point (TCP) [6], the critical chiral phase transition on the T-fi plane have 
been also attracting much interest. 

Recently, together with the energetic progress for the heavy-ion-collision (HIC) experiment facilities, such as the 
relativistic heavy-ion collider (RHIC) and large hadron collider (LHC), one can now probe hot and dense QCD matters, 
i.e. quark-gluon plasma (QGP), experimentally. It has been reported that very strong magnetic field in the order of 
the several times of to 2 [GeV 2 ] can be produced in the noncentral (peripheral) heavy-ion collision (HIC) experiment 
by STAR collaboration at RHIC [7, 8]. According to this strong magnetic field and CP- violating domains created 
inside the QGP, signals for possible P and CP violations were observed as the charge separation along the direction of 
the external magnetic field, which is perpendicular to the collision plane. Theoretically, this phenomenon is nothing 
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but the axial anomaly of electromagnetic currents [7-16]. The charge separation at relatively low T has been already 
investigated within the instanton- vacuum framework by one of the authors (S.i.N.) [17-19], in which related works 
and references can be found. Actually even before this sort of studies receiving much more interest recently due to 
the energetic progress of the HIC physics, QCD under magnetic field had been an important subject [20-24]. Inside 
QCD matter in the presence of the external magnetic field, the spins of the quarks are aligned along the direction of 
induced magnetic field according to their helicities. As a result, the quark-antiquark pair couples strongly, i.e. which 
is a phenomena denoted by the magnetic catalysis [21, 25]. Hence, taking into account that the order parameter of the 
spontaneous breakdown of chiral symmetry (SBxS) is the chiral (quark) condensate (qq), one can observe that SB\S 
will be enhanced by the presence of the external magnetic field. Accordingly there appear specific consequences as 1) 
the enhancement of the critical T for SBxS, T c , 2) the increase of constituent-quark mass, and 3) modifications in the 
low-energy constants (LEC). In the present work, we will address to all of these interesting consequences. To date, 
there have been many related works, for instance, from the Nambu-Jona-Lasinio (NJL) model [20, 21], holographic 
QCD [26, 27], lattice QCD simulations [28], linear sigma model [29], Polyakov-loop inspired model [30], and so on. 

In the present work, we want to investigate the (partial) chiral restoration under the strong external magnetic 
field in QCD matter. For this purpose, we employ the instanton-liquid model, modified by the Harrington-Shepard 
caloron solution finite T [31]. Although this approach does not manifest the quark confinement, such as the nontrivial 
holonomy caloron, i.e. the Kraan-van Baal-Lee-Lu caloron [32, 33], as shown in Refs. [17-19, 24, 34], it is a useful 
nonperturbatve method to study QCD matter at finite T. To include the external magnetic field, we make use of the 
linear Schwingcr method [17, 18, 35]. Besides we also take into account the meson-loop corrections (MLC) for the 
SU(2) light-flavor sector as the large-iV c corrections. MLC is essential to reproduce the correct current-quark mass 
dependence of relevant physical quantities [36] and universal-class chiral restoration pattern at finite T [24] . 

Our numerical results show that the chiral order parameters, such as the constituent-quark mass and chiral con- 
densate, are enhanced with respect to Bo due to the magnetic catalysis effect. There is a region where the u- and 
d-quark constituent masses coincide at e£?o ss (7 ~ 9)mJ, even for the explicit isospin-symmetry breaking, i.e. 
m u 7^ irid- The critical T for the chiral restoration, T c tends to be shifted higher pronouncedly in the presence 
of the B in the chiral limit. On the contrary T c keeps almost stationary for the physical quark mass case. The 
strength of the isospin breaking between the u and d quark condensates is also explored in detail by defining the ratio 
1Z = ((iu^u) — (iS d)) I ((iu^ u) + (id^d)) as a function of T and Bo- Finally, we compute the pion weak-decay constant 
and pion mass below T c as functions of T and B , showing correct partial chiral restoration behaviors. Our 
result also shows the changes of the F v and m„ due to the magnetic field, are relatively small in comparison to those 
caused by the finite T effect. 

The present work is structured as follows: In Section II, we make a brief introduction of the basic instanton-liquid 
model at vacuum without MLC and explain typical procedures to compute relevant quantities for further discussions. 
In Section III, we consider the inclusion of the MLC and the external magnetic field. The T-dependent modification 
of the instanton parameters using the Harrington-Shepard caloron is given in Section IV. Taking into account all the 
ingredients in the previous Sections, we derive the expressions for the saddle-point equation and chiral condensate as 
functions of T ad B in Section V. Section VI is devoted to presenting numerical results and associated discussions 
including the analysis of the pion properties at finite T and B . Summary and conclusion are given in Section VII. 



II. Ex A IN THE LEADING ORDER OF THE LARGE N c EXPANSION IN VACUUM 

We start by making a brief introduction for the instanton-liquid model in vacuum. This theoretical framework is 
characterized by the average of the inter-(anti)instanton distance R « 1 and that of the (anti)instanton size p w 1/3 
fm [37, 38]. The effective chiral action (Ex A) in the leading order (LO) of the 1/N C expansion can be written in 
Euclidean space as follows: 

Seff[m/] = C+N\nX + 2a 2 - J ^-^Tr^ln 

where C, J\f, A, <j, and m/ correspond to an irrelevant constant for further discussions, the instanton number density 
(packing fraction), the Lagrangian multiplier to exponentiate the 2Nf-t Hooft interaction, the saddle-point value of 
the chiral condensate, and current-quark mass for the flavor / for the SU(2) light-flavor sector. The Tr c j i7 indicates the 
trace over the color, flavor, and Lorentz indices. Detailed explanations on these instanton-related quantities are given 
in Refs. [37, 38]. In this picture quarks are moving inside the (anti) instanton ensemble and flipping their helicities. 
It results in that (anti)quarks acquire the momentum-dependent effective masses dynamically, i.e. constituent-quark 
masses. Assuming that the zero modes dominate the low-energy phenomena, we can write the Dirac equation for a 
quark for the (anti)instanton background as follows: 



ft + i[m f + M(k)] 
t + irrif 



(1) 



(ip + 4n) *jf = o, 



(2) 



3 



where Ajj and $// denote respectively the singular-gauge (anti)instanton solution and eigen function of the Dirac 
equation in the coordinate space [38] . By performing Fourier transformation of the <J? //, one is led to a momentum- 
dependent effective quark mass: 



M(fe) = M a = M f F 2 (k), F(k) = 2t 



/ (t)A-i(t) - h(t)K (t) - -I^K^t) 



t 



\k\p 



(3) 



where Mf stands for the constituent-quark mass for each flavor /, K n and I n are the modified Bessel functions [37]. 
Note that the F(k) can be interpreted as a quark distribution and plays the role of a natural UV regulator. Hence, in 
the instanton approach, UV divergences are regularized by construction without inserting any artificial form factors 
by hand. In practice, it is much easier to employ a parametrized form of the F(k) as in Refs. [17-19]: 



F(k) = 



2 + k 2 p 2 



(4) 



From the E%A we derive the following self-consistent (saddle-point) equations. They are used to determine relevant 
quantities such as the constituent-quark mass Mf at zero-momentum transfer k 2 = in Eq. (3): 



dS cS [m f ] _ dS cS [m f ] _ 
d\ ' da 

Similarly the chiral condensate can be computed by differentiating the E%A with respect to mf. 

1 dS eS [mf] 



Nf dmf 



(iqU)- 



(5) 



(6) 



From the first equation in Eq. (5), and Eq. (6), one obtains the expressions for the LO contributions for the instanton 
number density and the chiral condensate as functions of to/ as follows: 



A/"lo = 2N C N 



d A k 

WT 4 



M a M a 
k 2 + M 2 



LO 



4N r 



d*k 

WT 4 



M a 



TO/ 



k 2 + M 2 k 2 + ■ 



(7) 



Here we define M a = m f+M a . The value of Mlo is determined from the parameter R and its phenomenological value is 
~ (200 MeV) 4 [37]. At the chiral limit to/ =0, one solves the first equation of Eq. (7) self-consistently. The value of Mf 
turns out to be about 325 MeV. It is well consistent with the constituent-quark mass employed in usual quark models, 
i.e. 3Mf w M nuc i con . Applying this value of Mf into the second equation, we have (iq^q) = —(qq) ~ (235 MeV) 3 
for the SU(2) light-flavor sector. Again, this value of the chiral condensate is well matched with phcnomenologically 
accepted ones. 

In contrast to these seemingly successful numerical results for the chiral limit, the LO results for the Mf with finite 
to/ is considerably deviated from the available LQCD data [17-19]. In Refs. [36, 39], it was suggested that the correct 
to/ dependence of Mf can only be achieved by the inclusion of the meson-loop corrections (MLC) which is related to 
the next-to-leading order (NLO) of the large-./V c corrections. It has been also discussed that this NLO contributions 
play a critical role to reproduce the appropriate universal-class pattern of the chiral restoration as a function of T [19]. 
Consequently we will discuss the inclusion of the MLC and the magnetic field in the next Section. 



III. Ex A WITH MLC AND B FIELD 



Here we use a standard functional method [36, 39] to tackle the MLC corresponding to the large-A^ corrections. 
Taking into account the mesonic fluctuations around their saddle-point values, one can write the E%A via a standard 
functional method as follows: 



S cS [mf] = C+N\n\ + 2a 2 



J (27T) 



Tr 



c,/,7 



In 



h + iM a 



f, a +im f 



It 



d 4 k L_J_ 
(2^) 4 \ 4a 2 



LO 
d A k 



(27T) 



-Tr 



c,/,7 



M a 



M h 



% + iM a h + iM b 



(8) 
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where r, = (1, 75, it, 1x75) relates to the fluctuations from the isoscalar-scalar, isoscalar-pseudoscalar, isovector-scalar, 
and isovector-pseudoscalar mesons with the Pauli matrix denoted by r. As for the saddle-point values, we integrated 
out all the meson contributions except for the scalar one which signals for the SBxS. For more details for Eq. (8) can 
be found in Refs. [36, 39]. The k a and k b denote k and k + q, respectively. 

To study the impact of the external EM field on the QCD matter we need embed the external EM field in the E%A. 
Since we are only interested in the external magnetic field, assuming that it is static and aligned along the z axis as 
B = B z, we can choose the EM field configuration as follows: 



-Bo B Q 



(9) 



This EM field configuration in Eq. (9) makes the field-strength tensor, satisfying B = Fi 2 . Employing the linear 
Schwinger method [17, 18, 35], one can write the E%A in Eq. (8) as a function of the EM field strength: 



5 cff [m / ,F M1/ ] 



C + M In A + 2a 2 F - J T^Tr Ci/)7 In 



fia + %M A 



2 ^ 

1=1 



d 4 q 



\nl 1 



1 

4^ 



d A k 



-Tr. 



$ a + imf 
M, 

c,/,7 



% + iM a ft b + iM b 



(10) 



Here, the subscript F denotes the quantities under the external EM field. In Eq. (10), we use the notation = 
kp + efAp. Here e/ = Q/e, here e is the electric charge of the proton. Note that we have assumed that the NLO part is 
not modified by the external EM field, since the EM contribution from the NLO is small in comparison to the leading 
one. In the presence of the external EM field, the quark propagator is modified and approximated as [17, 18, 40]: 



1 _ p a -i[M a + e f N a (<r-F)] 
fi a + iM a ~ "~~ kl + M* 



0(Q"- 3 ) 



(11) 



where we have ignored the terms proportional to 0(Q^~ ), taking into account that (Q u ,Qd) = (2/3,-1/3). The 
relevant mass functions are also defined by the following expressions: 



M a =m f +M a = to/ + M f 



2 + fc2p2 



N a 



AM f p 2 



(2 + Agp 2 ) 3 " 



(12) 



From Eq. (10), the LO contributions for the instanton number density N in Eq. (5) and chiral condensate in Eq. (6), 
in the presence of the external magnetic field, can be derived up to 0(Q 2 ): 



LO,F 



2N c N f 



d 4 k 

WT 4 



+ 



2KB) 



(iq' ! q)hO,F = 4iV c 



d 4 k 



kl + Ml kl 



(13) 



where we assigned Bf as efBg = Qf(eBg) for simplicity, and the terms proportional to O(Qf) does not appear in the 
A/lo,f, due to Tr 7 (er ■ F) = 0. If B a = the expressions for the Af and chiral condensate in Eq. (13) recover those for 
vacuum given in Section II. Note that, although we do not have the explicit terms oc B 2 for the chiral condensate as 
far as we employ the quark propagator in Eq. (11), the condensate depends on the magnetic field because the Mf itself 
does implicitly. Moreover, the constituent-quark mass under the external magnetic field, which is determined from 
the saddle-point equation in the first line of Eq. (13), becomes different for the two flavors u and d, i.e. M u ,p ^ M^^f- 
It is because that they behave distinctively according to their electric charges. Similarly the chiral condensate also 
becomes flavor-dependent as shown in the second line of Eq. (13). Here is one caveat: For all the ingredients discussed 
so far, we have assumed that the instanton-packing fraction, M is immune from the external magnetic field as well as 
the flavor degrees of freedom by considering that the (anti)instanton is electrically neutral and non-flavored object. 

It is worth mentioning the differences between our theoretical framework and other chiral models. For instance, 
using the NJL model, Refs. [20, 21] obtained the magnetic-field dependent effective action. Most apparent difference 
between the present approach and theirs is how to regularize the UV divergence appearing in relevant physical 
quantities. The UV divergence is regularized by the nonlocal quark-instanton interaction in the present approach (see 
the quark form factors given in Eqs. (3) and (4)), whereas the regularization is achieved by adding and subtracting 
the lowest-Landau level (LLL) and the vacuum contribution of the chiral order parameters in the NJL model [20, 21]. 
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In this sense, in our approach, all the Landau levels are taken into account by construction in principle. Hence we 
do not need any specific analytic manipulations. A more detailed discussions on this regularization process based on 
the Landau level for the chiral models is given in Ref. [41] where the highest levels are naturally cut off, since the 
constituent mass at large momenta is suppressed when the magnetic becomes strong. 

The numerical value for of in Eq. (13) is obtained from the relation o 2 = Af/2 in the LO contributions as in 
Ref. [34]. Thus, using the LO part in the right-hand side of Eq. (13), we can write as follows: 

2 ■A/lO.F i att +- • *i ,. _ ■A/lo.f /,,n 

dp = — — h NLO contributions ss — ^— , (14) 

where we have rather safely ignored the NLO ones for the numerical calculations, since the NLO contributions are 
finite but much small in comparison to the LO one. By doing this, one can express the op simply as a function of the 
external magnetic field. Considering all the ingredients discussed so far, we write the J\f containing both of the LO 
and NLO (MLC) contributions, using Eq. (14): 

where the reduced functions Fi^ 3 are defined as 

M a M a 2N 2 M a M b [k a -k b + M a M b + M a M b + ^(M a + M b )} 
Fi{k) = — p-, i^ 2 (fc) = — , F 3 {k,q) = . (16) 

a a a, b 

Here, we have used the notations D\ b = k\ b + M 2 b . Similarly, the chiral condensate can be evaluated from Eq. (6) 
as follows: 

Here, the functions G\^z are assigned as 

n (u\ M a n ... m f M a M b {M a + M b ) 

G ^ k ) = "02 > G 2 (fc) = G 3 (fc,<z) jy^jyi , (18) 

a a b 

where D 2 = k 2 + m 2 . It is interesting to see from Eqs. (13) and (17) that, if the isospin symmetry is almost intact, 
m u ~ m di the difference between the condensates of the two flavors, i.e {iu^u) — (id^d), becomes negligible for the case 
with B = 0. However, as the strength of the magnetic field increases, the difference is also proportional to (e^ — e^)B 2 . 
For a better look on the isospin breaking effect, we define a quantity indicating the strength of the isospin breaking 
effect as follows: 

_ (ivtu) - (idU) 

~ (iutu) + (idU) ' { ' 

We also note that the ratio TZ is deeply related to the low-energy constant of the x?T Lagrangian, /i 3 [42, 43]. 



IV. INSTANTON PARAMETERS AT FINITE T 



To investigate the physical quantities in hand at finite T, we want to discuss briefly how to modify the instanton pa- 
rameters, p and R at finite T. We will follow our previous work [34] and Refs. [31, 44] to this end. Usually, there are two 
different instanton configurations at finite T, being periodic in Euclidean time, with trivial and nontrivial holonomies. 
They are called the Harrington-Shepard [31] and Kraan-van Baal-Lee-Lu calorons [32, 33], respectively. The non- 
trivial holonomy can be identified as the Polyakov line as an order parameter for the confinement-deconfinement 
transition of QCD. However, since we are not interested in the confinement-deconfinement transition in the present 
work, we choose the Harrington-Shepard caloron for the parameter modifications at finite T. We write the instanton 
distribution function at finite T with the Harrington-Shepard caloron as follows: 

d(p, T) = C Nc A^ s p N * p"- 5 exp [-(A N T 2 + hMp 2 ) P 2 } . (20) 

v ' 

c 
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Here, the abbreviated notations are also given as: 

0= -61n[A RS p cut ], p = -61n[A RS (i?>], C Nc 



4 60 e _1 - 68a Rs^Vc 
n*(N c -2)\(N c -l)\' 



Am = 



1 



11 



-Nr-1 



7 = 



27 



N? 



N 



(21) 



(22) 



Note that we defined the one-loop inverse charge $ and at a certain phenomenological cutoff value p cut and (i?) w R. 
As will be shown, only j3 is relevant in the following discussions and will be fixed self-consistently within the present 
framework. Ars stands for a scale depending on a rcnormalization scheme, whereas V3 stands for the three-dimensional 
volume. Using the instanton distribution function in Eq. (20), we can compute the average value of the instanton 
size, p 2 straightforwardly as follows [45]: 



p 2 (T) 



J d PP 2 d(p, T) _ [A%T* + AvfrAf] 5 - A N T 2 



Jdpd{p,T) ' IfoAf 
where v = (b — 4)/2. Substituting Eq. (23) into Eq. (20), the distribution function can be evaluated further as: 

1 



d(p, T)=C p"- 5 cxp [-M(T)p 2 ] , M(T) = -A N T 2 



(23) 



(24) 



The instanton-number density Af can be computed self-consistently as a function of T, using the following equation: 

Ar^M(T) = [CT(v)]i , (25) 

where we have replaced NT/V 3 -4- Af, and T(v) indicates a T function with an argument v. Note that C and can 
be determined easily using Eqs. (23) and (25), incorporating the vacuum values of the Af and p: C s=s 9.81 x 10~ 4 and 
j3 w 9.19. At the same time, using these results, we can obtain the average instanton size p as a function of T with 
Eq. (23). 

Finally, in order to estimate the T dependence of the constituent-quark mass Mf, it is necessary to consider the 
normalized distribution function, defined as follows: 



d N {p,T) - 



d(p,T) 
Jdpd(p,T) 



"'- 6 M V (T) exp [-M{T)p 2 



(26) 



Now, we want to employ the large- N c limit to simplify the expression of dM{p,T). Since the parameter b is in the 
order of 0(N C ) as shown in Eq. (21), it becomes infinity as N c -4- 00, and the same is true for v. In this limit, as 
understood from Eq. (26), dp{(p,T) can be approximated as a 5 function [37]: 



lim d N (p,T) = 5{p- p). 

N c —>cc 

Remember that the constituent-quark mass can be represented by [37] 

Mf oc VJ7 J dpp 2 S(p -p) = VAf p 2 , 
where Af and p are functions of T implicitly. We can modify Mf as a function of T as follows: 

f VA7 p 2 ^ 



Mf 



M f 



M f (T) 



(27) 



(28) 



(29) 



where Afo and po are those at T = 0. The numerical results for the normalized p/po and Af/Afo as functions of T arc 
in the left panel of Figure 1. As shown there, these quantities are decreasing with respect to T as expected: decreasing 
instanton effect. However, even beyond T* « Aqcd ~ 200 MeV, the instanton contribution remains finite. In the 
right panel of figure, we draw the quark mass as a function of T and absolute value of three momentum of a quark 
Ifel: 



M(k 2 ,T) = M f (T) 



2 + p 2 (T)k 2 



(30) 
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FIG. 1: (Color online) Normalized p/p and JV /JV as a function of T for iV c = 3, where M = N/V (left). The M(k 2 ,T) in 
Eq. (30) as a function of T and absolute value of the momentum |fe| (right). 



Note that we have ignored the Euclidean-time component of the four momentum by setting — 0. This tricky 
treatment simplifies the calculations in hand to a large extent, and we also verified that only a small deviation 
appears in comparison to the full calculations. Moreover, p in Eq. (30) is now a function of T as demonstrated by 
Eqs. (23) and (29) previously. As shown in the figure, M(\k\,T) is a smoothly decreasing function of T and |fe|, 
indicating that the effect of the instanton is diminished. Here, we choose Mj = 325 MeV at T = in drawing the 
curve as a trial. For more details, one can refer to the previous work [34]. 



V. Ex A WITH MLC AND EXTERNAL MAGNETIC FIELD AT FINITE T 

In this section we generalize our model to the finite temperature to explore the chiral restoration at finite T. For 
this purpose, we employ the Matsubara formula for the fermions. In this Euclidean-time description, it can be done 
by replacing the four-dimensional integral measure in the E^A into a three-dimensional one with a summation over 
the fermionic Matsubara frequency, w m = (2m + l)7rT: 

4>k Tff^L. ,31) 



(2tt) 4 ^ / (2tt 

in— —go 

Using this simple replacement we can rewrite the saddle-point equation in Eq. (15) as follows: 

N a 2N c N f f ^[TkW + J*(fc)B?] + P ^ ^ l I ' Jl (32) 

Summing over the Matsubara frequency, we define functions Ti as follows: 

T l =T s £ J F l - (33) 

m 

The analytic expressions for those functions are given in the Appendix. Note that the above equation recovers the 
expression of Eq. (16) in Ref. [24] when B = 0. Before going further, it is worth mentioning several assumptions made 
for deriving Eq. (32): 

• The instanton packing fraction, Af is a decreasing function of T indicated by the previous Section. However, 
we assume that the is not affected by the external magnetic field, since the (anti)instantons are electrically 
neutral. Moreover, the immunity of Af to flavors are assumed. 
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• In the momentum-dependent quark mass M a ^, we replace k 2 into fc 2 , ignoring the temporal term oc fcj = 
(2m + 1)ttT. We verified that this treatment makes the numerical calculations much simpler, and only small 
deviation was observed in comparison to the full calculations. Hence, the mass-related functions in Eqs. (4) and 
(12) are redefined as follows: 

AM f 4M f p 2 
a ~ (2 + fcV) 2 ' a_ ~(2 + fcV) 3 ' ( ' 

As for the M b and N b , we replace k with k + q. 

• Similarly, the term k a ■ kt = k 2 + k ■ q in the functions of Fi and is replaced by + k 2 + k ■ q. Moreover, 

+ 



the denominator is also replaced by D 2 b = + k 2 b + M 2 



• We also replace the integral variable (74, which corresponds to the fourth-component of the pion momentum, 
into an additional parameter A as in Ref. [36]. Since the isovector-pseudoscalar meson, i.e. pion dominates the 
meson fluctuations, it is reasonable to set the cutoff A proportional to m T as follows: 

A^m^. (35) 
P 

Note that, in the above equation, we have multiplied a factor po/p to m^ in order to include T dependence of 
the cutoff mass. Moreover, this multiplication factor represents a correct chiral-restoration pattern of m^, i.e. 
the mass of the pion, as a Nambu-Goldstone (NG) boson, increases as SB^S restored partially. 



Similarly the chiral condensate in Eq. (17) reads: 
Again, the relevant functions Qi = TJ2 m @i are defined and given in the Appendix. 



VI. NUMERICAL RESULTS 



We present and discuss our numerical results in this Section. We choose R w 1.0 fm and p w 0.34 fm which give 
Mf « 315 MeV in the present framework [43]. The values for the current-quark mass are reported as m u = (1.7 ~ 3.3) 
MeV and nid = (4.1 ~ 5.8) MeV [46]. Thus, we take the average values, (m u ,md) ~ (2.5,5) MeV. For simplicity 
we introduce an positive integer n. The magnetic field is assigned in terms of the pion mass as ei?o = nm 2 ,. In the 
Gauss unit for the magnetic field, we have the convention, B ~ n (1.2 x 10 18 ) G. For instance, n w 1 corresponds to 
a neutron star or magnetar with very strong magnetic field. The case with n w 10 or more may be observed inside 
the quark-gluon plasm created in the ultra high-energy peripheral heavy-ion collisions, such as RHIC and LHC, as 
the main source for the nontrivial QCD vacuum effect, i.e. chiral magnetic effect [7-16]. 



A. Constituent-quark mass for each flavor: Mf 

Here we present the numerical result of the T and B dependencies of the constituent- quark mass Mf which is 
one of the order parameters for the (partial) chiral restoration. According to the universal class of the restoration 
pattern, as for the chiral limit m Uy d = 0, the restoration pattern shows the second order, whereas it becomes the 
crossover for the case with the physical quark mass, i.e. m u ^d = (2.5,5) MeV. In our previous work [24], the MLC 
contributions, as the large- N c corrections, play a critical role to reproduce the universality in an appropriate manner. 
Similar observation was also reported in Ref. [47]. Here we study the impact of the external magnetic field on the 
two light-flavor QCD matter at finite T. 

In Figure 2, we depict the results of the Mf at T = for the chiral limit in the panel (A) and physical quark mass 
case in the panel (B). In the absence of the magnetic field (n = 0), we have M Ut< j, = 315.02 MeV for the chiral limit, and 
M u _d = (316.16,317.04) MeV for the physical quark mass. Accounting for m u < raid and that the constituent-quark 
mass becomes M u> d ~ Mf + rrif, the observed results can be easily understood. As the external magnetic field emerges 
the constituent-quark mass becomes heavier. It is due to the magnetic catalysis. Since the effects of the magnetic 
catalysis is proportional to e 2 as in Eqs. (32) and (36), the M u grows rapidly more than the Md with respect to the 
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external magnetic field, i.e. > e^. Beyond n = (17 ~ 18), the curve for the M u starts decreasing slightly. It is 
interesting to see that, for the physical quark mass case in the panel (B), there appears a point at which the M u and 
M d coincide each other (n f«7.5). It is because that the broken isospin symmetry (m u ^ m d ) is compensated by the 
effect from the external magnetic field. From our results, this interesting phenomena appears at the magnetic field 
Bo w 10 19 G, which can be created at the heavy-ion collision experiments. Beyond this point, the ordering of the u- 
and (i-quark constituent masses are reversed. 

Here, we want to discuss briefly the saturation behavior of the constituent-quark mass observed in the panel (A) 
and (B) beyond n w 15 for the M u . This can be understood as follows: As discussed in the previous Sections (see 
Eq. (11) for instance), we have taken into account the contributions up to 0{Q 2 ). Hence, as for the stronger magnetic 
fields, higher-order contributions which are not included in the present work will be no longer negligible. Accordingly, 
we verified that this saturation behavior becomes weaker in the presence of higher-order contributions. Nevertheless, 
the complete and systematic treatment of those higher-order terms is left for the future work. 

In the panel (C) and (D) the results for the case at T = 50 MeV are demonstrated. Note that, for all the cases, 
the absolute values for the Mf decrease in comparison to those at T = 0. But the shapes and behaviors of the curves 
are similar. This decreasing tendency can be understood by the diluting instanton ensemble at finite T. For more 
details on the diluting ensemble in the present framework, one may refer to Refs. [34]. It turns out that the decreasing 
rate of the constituent-quark mass is a few percent from T = to T = 50 MeV: M u ,d = (302.96) MeV for the chiral 
limit and M u ^ = (304.44, 305.62) MeV for physical quark mass, at n = 0. Similarly to the vacuum case T = 0, 
there appears a point n w 8, at which the u- and d-quark constituent masses are very close to each other. Moreover, 
the isospin breaking effect (m<j ^ rn u ) becomes more obvious at T = 50 MeV compared to that for T = 0. It is 
because of the decreasing dynamically-generated quark mass at finite T. All of these observations are all based on 
a nontrivial competition between the magnetic catalysis and diluting instanton effects at finite T: The former tends 
to enhance the constituent quark mass, whereas the latter suppress it. Hence, the position of the equal mass point 
is a consequence of this nontrivial competition between the two mechanisms, on top of the explicit isospin symmetry 
breaking. The reversing order of mass beyond the equal mass point is also observed at finite T. In Ref. [21], the 
authors found qualitatively the same results; monotonically increasing curves for the Mf with respect to the external 
magnetic field was observed and the M u is more sensitive to it. However, the rate of increasing is much higher than 
ours. 



B. Chiral condensate: (iq'q) 

Now we are in a position to discuss the (partial) chiral restoration in the presence of the external magnetic field. 
It is indicated by the chiral order parameter, i.e. quark condensate. In Figure 3, we show the numerical results of 
the chiral condensate for the d (left column) as well as u (right column) quarks, separately. We observe the second 
order chiral phase transition for the both flavors in the chiral limit case shown in the panel (A) and (B). This result is 
expected from the universal class of the restoration pattern. Note that this correct restoration pattern in the present 
instanton framework is only achieved by the inclusion of the MLC as the large- N c corrections [24]. Turning on the 
external magnetic field, one finds that the SB%S is enhanced, that is, the values for the quark condensate and T c both 
increase for the two flavor. Among them we see that the u quark condensate is more sensitive with respect to the 
magnetic field. It is due to the larger quark electric charge of the u quark. The critical T for the both flavors, T" 
and T c d are listed in Table I. We observed that {iv)u) w (icftd) w (247 MeV) 3 at T = 0, which is just compatible to 
the empirical value of the isospin-symmetric quark condensate about (250 MeV) 3 . 

Considering the physical quark mass case, the chiral phase transition for the two flavors are shown in the panel 
(C) and (D) of Figure 3. Following the universal class of the restoration pattern, the curves represent the crossover. 
The magnetic field effects are negligible for the d-quark condensate in the panel (C), due to the smaller quark electric 
charge, whereas w-quark condensate in the panel (D) shows visible changes in the vicinity of T w 180 MeV with 
respect to the magnetic field. The T c can be obtained by computing the inflection point of the curves for the crossover 
phase transition [48], resulting in that T c d w 200 MeV for all the n values and T™ = (180 - 200) MeV for n = (0 ~ 20). 
It is worthy of noting that the changes in the T c for the physical quark mass, due to the magnetic field, are relatively 
small in comparison to those for the chiral limit. This tendency is qualitatively consistent with the lattice QCD 
estimations [49]. As for the physical quark mass case, the quark-condensate values for the both flavors at T = are 
about (248 MeV) 3 , which almost coincides with those for the chiral limit. 
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FIG. 2: Constituent-quark mass (Mf) as a function of n = eBo/rri^ for the chiral limit (left column) and physical quark mass 
(right column). In the upper and lower rows, we show the numerical results for T = and T = 50 MeV, respectively. For more 
details, see the text. 





n = 


n = 10 


n = 20 


rpU 
C 


170.4 MeV 


173.9 MeV 


183.1 MeV 


rpd 

■*■ c 


170.4 MeV 


171.3 MeV 


174.1 MeV 



TABLE I: Critical temperature for the u and d flavors, T" and T%, for the chiral limit in the presence of the external magnetic 
field, n = eBo/ni^- 



C. Ratio of the two-flavor quark condensates: 1Z 



In the previous Subsections, we discussed the competition between the magnetic catalysis and diluting instanton 
effect at finite T, on top of the explicit isospin breaking. In the present Subsection, we want to take a more careful 
look on the isospin breaking of the quark condensates by defining a quantity as in Eq. (19). We also note that the 
ratio 1Z is deeply related to the low-energy constant of the xPT Lagrangian, /13 [42, 43]. In Figure 4, we depict the 
1Z for the chiral limit in the panel (A) and physical quark mass case in the panel (B). In the chiral limit without the 
external magnetic field, the u- and d-quark condensates are the same so that 1Z = for any T values. As the magnetic 
field increases from n = 0, the 1Z becomes a positive and stiffly increasing function. It is because that the u-quark 
condensate increases more rapidly than that for the d quark and the magnetic catalysis effect is proportional to as 
in Eq. (36). At the critical T, the values for the 1Z diverges, signaling the second-order chiral phase transition. 
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FIG. 3: Chiral condensate as a function of T. We draw the numerical results for the chiral limit and physical quark mass in 
the left and right columns. The results for the d and u quarks are given in the upper and lower rows, respectively. For more 
details, see the text. 



The situation becomes quite different for the physical quark mass case in the panel (B). Without the magnetic field 
the TZ decreases then becomes negative beyond T sa 100 MeV. As the magnetic field increases, the curves are shifted 
to higher T, and there appears a bump around T — 180 MeV for n = 20. The apparent difference between the two 
cases is not hard to understand. In the chiral limit, there is as a sort of intact degeneracy between u and (i-quark 
condensates. Such that nonzero TZ values are only possible in the presence of the finite external magnetic field which 
breaks this degeneracy [50]. This is also true for the nonzero degenerated quark mass of the flavors, m u = =/= 0. 
On the contrary, if this degeneracy between the quark condensates is lifted up beyond T ss 100 MeV with decreasing 
nonperturbative effect (instanton), the explicit isospin breaking becomes more pronounced and results in the negative 
difference of the condensates at the zero magnetic field. It is due to the fact that the heavier quark causes the larger 
quark condensate in general. On the hand the u quark is more sensitive to the magnetic catalysis because of its larger 
electric charge. In other words, the explicit isospin breaking effect pushes the TZ downward but the magnetic catalysis 
pushes the TZ upward. This competition also causes a bump around T — 180 MeV for the 1Z. Nevertheless TZ goes 
down when T increases beyond T = 180 MeV indicating the explicit isospin breaking effect due to the quark mass 
difference wins over magnetic catalysis there. 
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FIG. 4: (m, d)-quark condensate ratio 1Z in Eq. (19) as function of T for different n = eBo/m\ values. We present the numerical 
results for the chiral limit and physical quark mass in the left and right panels, respectively. 



D. Pion properties at finite T under the magnetic field: and 

Finally we want to make an analysis on the pion properties such as the pion weak-decay constant and pion mass, 
at finite T in the presence of the external magnetic field, below the critical T. Since in the Nambu-Goldstone phase 
(T is below 100 MeV), the isospin symmetry is only slight broken as shown in Figure 4. Hence, in what follows, we 
focus on the properties with the isospin symmetry. For this purpose we employ the Gell-Mann-Oakes-Renner (GOR) 
relation defined as: 

TO - = Yl J&fafQf) -* ^■<*« t ?)m / =o ■ (37) 



/= 



isospin symmetric 



Here, we defined fhf — (m u + m^)j1. Since the quark condensates have been already computed as functions of 
T as well as Bq in the previous Sections, it is enough to calculate the pion-weak decay constant in the same 
framework. For simplicity, we ignore the MLC contribution to compute the F n and the nonlocal contribution for the 
time being [51]. Then the analytical expression for the F v reads in the instanton framework for the vacuum as follows: 



d 4 k M%-k 2 M k N k 
{2ttY (fc 2 + M 2 ^ 2 



K « ^N c j— A * ■ J— , (38) 



where we again have assumed the isospin symmetry. The r\ denotes a correction factor for the case without the 
nonlocal contribution. From Ref. [51], the value for the r\ can be estimated as about 0.5 to obtain the empirical 
value for the pion-weak decay constant, i.e. F v ss 93 MeV. Although some dynamical information from the nonlocal 
contributions are missing by this simplification, it is still useful for a simple and qualitative analysis. If we induce 
the EM field externally, we can replace the constituent-quark mass squared approximately as M\ — > M| + 2 V|Z? 2 , 
according to Eq. (13). Moreover, taking into account that the term M^Nk can be obtained by differentiating Mj^/2, 
we have Mj^^k ~ M^N^ + 2N k (dNk/dk)B'j. Hence, we can write the expression for the F n as a function of T and 
-Bp, employing the fermionic Matsubara formula: 



- w E J ^3 (< + fc2+M|)2 - 4 ^ J ^ + ** + ^ ■ ( 39 ) 

Here, we have defined a flavor-averaged external magnetic field, i.e. the £? 2 = \ [B\ + Bj) , considering the isospin 
symmetric matter. Analytic expressions for the relevant functions JC\^3 are given in the Appendix. 

In Figure 5, we preset the numerical results of the pion weak-decay constant F w (A) and pion mass (B) as 
functions of T and the strength of the magnetic field. In our numerical calculations, we have chosen 2m j « 10 MeV 
in Eq. (37) as a trial, although we evaluated the analytical expression for the F„ near the chiral limit. In the panel 
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FIG. 5: Pion weak-decay constant F n (A) and pion mass m w (B) as a function of T, varying the strength of the magnetic field. 

(A), the F n smoothly decreases with respect to T indicating the partial chiral restoration. At T = 100 MeV, the 
value of the is about 10% reduced. Increasing the strength of the magnetic field one finds that the value of is 
enhanced by a few percent. At T = 0, we have F^ = (93.47,93.81,94.16) MeV for n = (0, 10,20), respectively. The 
effect from the magnetic catalysis appears more important in the higher T region. 

The numerical results for the m w are given in the panel (B). As a signal for the partial chiral restoration, the pion 
mass increases with respect to T but decreases with respect to Bq. In other words, the enhancement of the SB^S due 
to the magnetic catalysis is quite small. Only about 0.5 MeV decease in the pion mass is observed for n = (0 — > 20). 
Again, the magnetic catalysis plays an more significant role in the higher T region. 

VII. SUMMARY AND CONCLUSION 

In the present work, we have investigated the (partial) chiral restoration at finite T in QCD matter for the SU(2) 
light-flavor sector, under the strong and static external magnetic field. To this end, we employed the T-modificd 
instanton-liquid model, together with the linear Schwinger method and fermionic Matsubara formula. We also took 
into the meson-loop corrections as the large- N c corrections to reproduce the correct phase transition pattern. We 
then present the numerical results for the constituent-quark mass, chiral condensate, and isospin-symmetry breaking 
effect as functions of T, B , and flavor degrees of freedom. Below, we list important theoretical observations of the 
present work: 

• Relevant instanton parameters R and p are modified as functions of T, resulting in the diluting instanton 
ensemble with respect to T, i.e. decreasing the SBxS effects. The external magnetic field enhances the SB%S in 
terms of the magnetic catalysis, which is proportional to (efB ) 2 . Hence, the u-quark constituent mass is more 
sensitive to the magnetic field and increases considerably more with respect to the field strength compared to 
the d quark. 

• On top of the explicit isospin symmetry breaking, there appears a point at which the constituent quark masses 
for the u and d quarks coincide each other for the strong magnetic field eBo 10 19 G. In the chiral limit, we 
observe the second-order chiral phase transition, as expected from the universal restoration pattern, effected 
much by the meson-loop corrections. Naturally, the crossover phase transition takes place for the physical quark 
masses. 

• The effects from the magnetic catalysis becomes more pronounced in the higher T region because as T increases 
the SB%S effects generated from the instanton is weakened so that the magnetic catalysis effects becomes 
relatively more important. The critical T, i.e. T c is shifted to higher T due to the magnetic catalysis, whereas 
the change of the chiral condensate values is relatively small. T c becomes flavor-dependent because the magnetic 
catalysis effect depends on the electric charge of the quarks. 
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• The isospin breaking between the quark condensates is explored by defining the ratio 71 as a function of T. 
As for m u = rrid case the ratio is zero at Bq = and monotonically increases with respect to T for the finite 
magnetic field. For the physical quark mass case, the ratio 1Z shows nontrivial structures with respect to T and 
B n due to the complicated competition between the magnetic catalysis and the explicitly isospin breaking effect 
which becomes more important at higher T because of the decreasing nonperturbative effects. 

• According to our simple and qualitative analysis using the GOR relation, we observe correct partial chiral- 
restoration and magnetic-catalysis behaviors for the pion-weak decay constant and pion mass m v . They 
decreases and increases about 10% atT« 100 McV in comparison to those at T = 0, respectively. However, 
the changes due to the magnetic field are relatively small, just a few percent. 

Now we obtain an effective chiral action at finite T as well as the magnetic field for the physical quark mass. 
If the strong magnetic field is created in the peripheral heavy-ion collisions as reported, it is worthy of studying 
the hadronization processes in the presence of the magnetic field, i.e. the dilepton production via the vector-meson 
dominance under the magnetic field for instance. Moreover, the QCD phase diagram on the fi-T plane and critical 
values, such as the critical end point (CEP) and tricritical point (TCP), are also able to be explored in our model in 
principle. Inclusion of finite \i to the effective action in the instanton framework is under progress and related works 
will appear elsewhere. 
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Appendix 

The relevant functions in Eq. (32) and (36) are given as follows: 

Tx = M a M a U 2 , T 2 = 2N 2 a H 2 , Fz = M a M b Hi + M a M b [k ■ (k + q)} H 3 , 



Fa 



M a M b 



M a M b + M a M b + ^(Ma + M b ) 



n 3 



Gi = M a n 2 , G 2 = -m f H u 5 3 = M a M b (M a M h )H 3 , 

K x = (M*-k 2 M a N a )H 5 , JC 2 = -M a N a 7i 6 , JC 3 = 2N 2 a % b . 



Hir^Q are explicitly given as follows, 
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= T E 
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w^ + El 2E a \2T 
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